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Abstract 

A self-adjoint operator with dimensions of time is explicitly constructed, and it 
is shown that its complete and orthonormal set of eigenstates can be used to define 
consistently a probability distribution of the time of arrival at a spatial point. 



I. INTRODUCTION 



The role that time plays in quantum mechanics has always been controversial. This 
is in part a consequence of the rather singular status that time exhibits in nonrelativistic 
physics. In particular, time enters the Schrodinger equation as an external parameter 
and accordingly the quantum formalism is usually concerned with probability distri- 
butions of measurable quantities at a definite instant of time. However, one may also 
ask for the instant of time at which a certain property of a quantum system takes a 
given value. In this case time has the character of a dynamical variable: It depends 
on the initial state of the system and on its dynamical evolution, and appears as an 
intrinsic property of the physical system under study. Since such an instant of time 
is, in principle, a perfectly measurable quantity it seems natural to try to incorporate 
the concept of a time observable into the quantum formalism. However this is not 
an easy task. The standard quantum formalism associates measurable quantities with 
self-adjoint operators acting on the Hilbert space of physical states, and postulates that 
the probability distribution of the outcomes of any well-designed measuring apparatus 
can be obtained in terms of the orthogonal spectral decomposition of the corresponding 
self-adjoint operator, with no explicit dependence on the particular properties of the 
measuring device. Therefore, the problem reduces, in principle, to find a suitable quan- 
tum operator. This is usually accomplished via the correspondence principle, starting 
from the corresponding classical expressions and quantizing by using certain specific 
quantization rules. However, in doing so one frequently has to face with the problem 



that in general there exists no unique way to obtain a quantum operator which reduces 
to a given known expression in the classical limit {h ^ 0). 

Given the Hamiltonian H{q,p) of a conservative classical system, expressed in terms 
of canonical variables {q,p), one can always make a canonical transformation to new 
canonical variables {H,T), where H is the Hamiltonian of the system and T its conju- 
gate variable, which satisfies the Hamilton's equation ||l],^ 

= (1) 

{H,T} denoting the Poisson bracket of H and T. The important point is that the 
above equation clearly reflects that the canonical variable T is nothing but the interval 
of time. Thus the next step would be to take advantage of this desirable fact and 
translate the above formulation to the quantum framework. This can be easily ac- 
complished by means of the canonical quantization method Q, which basically states 
that the classical formulation remains formally valid in the quantum domain with the 
substitution of Poisson brackets by commutators, {H,T} — > 1/ih [H,T], and interpret- 
ing the dynamical variables as self-adjoint operators in the Heisenberg picture. Then, 
based on the correspondence principle and the canonical quantization method one is 
led to look for a self-adjoint time operator conjugate to the Hamiltonian, 

[H,f]=ih. (2) 

This commutation relation, which as can be easily verified also holds true in the 
Schrodinger picture, has the additional desirable consequence that it implies the un- 
certainty relation 

AH AT > 1/2 \[H,f]\, (3) 

with AH and AT being the usual root-mean-square deviations of the corresponding 
dynamical variables. Unfortunately no such a time operator exists. As remarked by 
Pauli the existence of a self-adjoint operator satisfying the above commutation relation 
is incompatible with the semibounded character of the Hamiltonian spectrum [Q]. 

The lack of a proper time observable has a number of consequences In particular 
the time-energy uncertainty relation has remained unclear over the time. This is so 
basically because contrary to what happens with the well-known position-momentum 
uncertainty relation, there exists no unique way to put in a quantitative setting what 
is really meant by the time spread AT. In fact the consequences derived from incorrect 
application of the time-energy uncertainty relation have led to a great deal of confusion. 

Another related problem which remains controversial at present is that concerned 
with the formal definition of traversal and tunneling times §,0]. This subtle question 
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has received considerable attention in recent years |8-15| motivated, in part by the 
possible applications of tunneling in semiconductor technology. However, the simplest 
problem involving time as a dynamical variable is that concerned with the time of arrival 
of a free particle at a given spatial point. Such a time constitutes a well-defined concept 
which at a classical level can be extracted from the formalism by simply inverting the 
corresponding equations of motion. Moreover it is a perfectly measurable quantity 
whose probability distribution can, in principle, be experimentally determined within 
any desirable precision. However, the standard quantum theory of measurement does 
not provide any formulation which allows to infer such a probability distribution. In 
fact, time ago Allcock [16| argued against such a possibility. This author claims that it 
is not possible to construct any operationally meaningful and apparatus-independent 
probability formula. Even though more recently a number of works addressing this 
question from a more optimistic perspective have appeared [pj.lT-p^, the problem is 



yet far from being resolved, and additional investigation on this fundamental question 
is worthwhile. 

In this paper we analyze the possibility of defining a probability distribution for 
the arrival time of a quantum particle at a definite spatial point. Specifically, we 
are interested in searching for an apparatus-independent theoretical prediction for the 
probability distribution of arrival times at a given spatial point, as a certain function 
of the initial state of the system. Our results turn out to be similar to those previously 
obtained by Kijowski |17|. However, the approach by Kijowski is based on the definition 
of a non-conventional wave function which evolves on a family of x = const planes 
(instead of evolving in time according to the Schrodinger equation) , and whose relation 
to the conventional wave function is unclear |2^,24|. Our approach, on the contrary, is 
entirely developed within the formalism of standard quantum mechanics. 

We begin considering in detail the case of a free quantum particle and then we 
study the more interesting case of a quantum particle scattered by a potential barrier. 
For simplicity we shall restrict ourselves to a one spatial dimension. For our purpose, 
it proves to be useful analyzing first in some detail the reason for the nonexistence of 
a self-adjoint time operator in quantum mechanics. 



II. NONEXISTENCE OF A TIME OPERATOR IN 
QUANTUM MECHANICS 

As stated above, according to Pauli's argument, because of the semibounded char- 
acter of the energy spectrum there exists no self-adjoint operator conjugate to the 
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Hamiltonian, i.e., satisfying the commutation relation @. The same negative con- 



clusion was found by Allcock |16| using a somewhat different argument based on the 
time-translation property of the arrival time concept. 

If {jr)} denotes a set of measurement eigenstates for the arrival time at a given 
spatial point of a particle in the quantum state then according to the standard 
quantum formalism the probability amplitude for the arrival time at the instant t = T, 
would be given by il^(T) = {Tltp). If one translates the state of the system forward 
through time by an amount r, i.e., 1-0) — > IV'') = exp{—iHT/h)\'tp), then it seems 
natural to expect the probability amplitude to transform according to ip{T) — > ip'{T) = 
ip{T + t). Since this transformation property must be true for any state vector \ip) it 
follows that the measurement eigenstates must satisfy 

\T + t) =e'^^/^\T), (4) 

which reflects the fact that under a translation backward in time by an amount r, any 
measurement eigenstate corresponding to arrival time at the instant t = T, transforms 
into another measurement eigenstate, corresponding to an arrival time t = T + t. 
Based on general grounds, Allcock showed that measurements eigenstates with such 
a desirable property cannot be orthogonal, which implies that it is not possible to 
construct the corresponding self-adjoint arrival-time operator. It is not difficult to see 
that this negative conclusion can be traced back again to the semi-infinite nature of 
the Hamiltonian spectrum. To this end let us consider the following three statements: 

i) There exists a self-adjoint operator T conjugate to the Hamiltonian i.e., sat- 
isfying [H,T] = ih. 

ii) There exists a self-adjoint operator T, whose (orthonormal and complete) set of 
eigenstates transforms under time-translations as e^^'^^^ \ T) = |T + r). 

iii) There exists a self-adjoint operator T which generates unitary energy-translations, 
i.e., such that for any energy eigenstate \E) and any parameter e with dimensions 
of energy, it holds 

e^^'/^\E) = \E-e), (5) 

where the operator T is assumed to be defined onto the whole Hilbert space 
spanned by the Hamiltonian eigenstates. 

It is not difficult to see that these statements are in fact equivalent. Indeed, if i) is 
true, then by induction, one has 

[H'',f]=inhH''-\ n>l, (6) 
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where = 1. Of course the vahdity of Eq. rests on the reasonable assumption 
that the Hamiltonian is weh-behaved enough as to guarantee the existence of all its 
higher integer powers. Since it also holds that [H^,T] = for n = 0, then multiplying 
Eq. @ by (ir/h)"' /n\ (r being an arbitrary parameter with dimensions of time) and 
summing from n = to n = oo, one finds 

[e'^^/^,f] = -Te'^^/^. (7) 

If {|r)} denotes a complete and orthonormal set of eigenstates of T, then according to 
Eq. (0) it holds 

fe'^^^^\T) = {T + T)e'^^^^\T), (8) 

which after suitable choice of normalization and phase leads to statement ii). Con- 
versely, if ii) is true for any eigenstate |r) and any parameter r, then one can repeat 
the same steps backward to reach i). 

On the other hand, it can be readily seen that statement i) also implies statement 
iii). Indeed, if i) holds one has by induction that 

[F,f"] = m^f"-\ n>0, (9) 
(r*^ = 1), which implies that for any parameter e with dimensions of energy 

[H,e''^'/^] = -ee'^'/^. (10) 
Therefore, according to Eq. (^) any energy eigenstate \E) verifies 

He'^'/^\E) = {E -e)e'^^/^\E), (11) 

from which after proper normalization follows iii). An analogous reasoning can be 
repeated from iii) to i), which shows the equivalence among the above three statements. 

Since iii) is obviously incompatible with a semibounded Hamiltonian spectrum, it 
follows that it is not possible to find a self-adjoint arrival time operator satisfying the 
desirable conditions i) or ii). 



III. A SELF-ADJOINT OPERATOR WITH 
DIMENSIONS OF TIME 

We start by considering the simplest conceivable arrival time problem, namely a 
one-dimensional free particle moving along the x-axis towards a detector. In looking 
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for a probability distribution of the time of arrival it is most convenient to work in the 
energy representation {\E,a); E >0, a = +, —}, defined by the eigenvalue equations 

Ho\E,±)=E\E,±), (12) 

P\E,±) =±V2mE\E,±), (13) 

where = P^/2m is the Hamiltonian of the free particle, and P its momentum 
operator. The orthonormal and complete set of energy eigenstates a)}, which 
satisfy 

/■oo 

J2 / dE\E,a){E,a\ = l, (14) 

{E,a\E',a')=Saa'S{E-E'), (15) 

can be expressed in terms of the usual momentum representation by means of the 
relation 

\E, ±) = {mjlE)^!^ \p = ±\/2mE), (16) 
where the momentum eigenstates {\ p)} are normalized as 

f+OO 



/" 

J — ( 



dp\p){p\=l, (17) 



{p\p')=S{p-p'). (18) 

As stated above, the impossibility of finding a time-of-arrival operator can always 
be traced back to the bounded character of the Hamiltonian spectrum. To circumvent 
such a difficulty we shall instead look for a self-adjoint operator T with dimensions of 
time, conjugate to a conveniently defined self-adjoint operator H, with dimensions of 
energy and a non-bounded spectrum, 

[H,f]=ih. (19) 

Of course this is a somewhat arbitrary procedure, since the definition of T depends 
in a fundamental way on the arbitrary choice one makes for the operator H. Moreover, 
as long as H differs from the Hamiltonian of the system, the corresponding opera- 
tor T could not be associated to the actual physical time. Therefore, it remains the 
fundamental question of verifying whether it is possible to give a proper physical inter- 
pretation to the selected T-operator in terms of measurements results, i.e., whether it 
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is possible to define an algorithm which enables us to connect the probability distribu- 
tion of measurement results with the set of eigenvalues and eigenstates of the operator 
T. At this point we shall postpone this essential question and simply consider the 
procedure just outlined to be worth exploring. 

We now introduce projectors, 0(±P), onto the subspaces generated by plane waves 
with positive/negative momenta, 

G(±P) = / dp\±p){±p\, (20) 
Jo 

and define the self-adjoint operator 

sgn(P) = G(P) - G(-P). (21) 
Obviously, sgn(P) commutes with the Hamiltonian and satisfies the eigenvalue equation 

sgn{P)\E,±)=±\E,±). (22) 

This operator allows us to define a simple self-adjoint operator with dimensions of 
energy, 

H = sga{P)Ho, (23) 
which exhibits a non-bounded spectrum, 

?t:|^,±) = ±^|^,±), (^>0). (24) 

Notice that this is, in a sense, the simplest choice, since the restrictions of Ti to 
the subspaces spanned by plane waves with positive/negative momentum coincide with 
plus/minus the corresponding restrictions of the Hamiltonian Hq. Specifically, 

G(±P) n Q(±P) = ± Q(±P) Ho Q(±P). (25) 

Introducing the following notation for the energy eigenstates 

i\+E)^\E,+) iie>0 
\^) - I \-E) = \E,-) if £<0 ^"^^^ 

the above results can be rewritten in terms of the complete and orthonormal set of 
states {|£); £ G (— cxd, -|-(X))} satisfying the eigenvalue equations 

H\e)=s\s), (27) 



P|£) = sgn(£)^2m|£| |£), (28) 
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Hole) = \e\\e). (29) 

Now searching for a self-adjoint operator T conjugate to is a straightforward 
matter. To this end let us introduce the states |t) defined in the { |e)} basis as 

/ + 00 
dee''^/^\e). (30) 

These states also constitute a complete and orthonormal set. Indeed, 

{t\t')= de{T\e){e\T') = h-U dee-''^^-^^'^ = 5{t - t'), (31) 

J—oo J—oo 



/+00 /■+00 r+oo r+co 

dT\T){T\ = h-^ dr d£ de' e'^^-'^^/^\e){e'\ 
~oo J—oo J~oo J — oo 

+ 00 

de\e){e\ = l. (32) 

We can therefore define a self-adjoint operator, with eigenstates and eigenvalues given 
by |r) and r, respectively, in terms of its spectral decomposition, 

/+00 
dTT\T){T\. (33) 
-oo 

The operator so defined has dimensions of time and automatically satisfies the com- 
mutation relation ([l9|) . However, there exists no guarantee that it will be useful in the 
time-of-arrival problem. In fact, T turns out to be invariant under time reversal, and 
consequently the variable r cannot be identified with the physical time t. This can be 
most easily seen in the momentum representation, 

{p\t) = (|p|/m/l)l/2 gisgn(p)i4r/ft_ (34) 

Let TZ denote the time-reversal operator, then we have 

/+00 r+oo 
dp \-p){p\t)* = dp \-p){-p\t) = |r), (35) 
-oo J —oo 

so that, according to Eq. (33), it holds that TZTW = T . Moreover, 

|r + r') = e*^"'/'* |r) / e''^°"'/^ |t), (36) 
and the states { |t)} do not exhibit the desirable time-translation property (H) either. 
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In spite of these facts, it is possible to give a physical interpretation to the operator 
T. As we shall see below, one can consistently define a probability distribution of 
arrival times in terms of the eigenvalues and eigenstates of this operator. To this end it 
is convenient to decompose |r) as a superposition of negative- and positive- momentum 
contributions, 

|r) = /i-^/M dEe-'^^/^\E,-) + h-^l'^ dE e'^^'^ \E,+) . (37) 
Jo Jo 

Defining new states 

roo 

\t,±) = h-^/^ dEe'^'/^\E,±), (38) 
Jo 

we see that |r) can be written in the form 

\t) = \t = -T,-) + \t = +T,+). (39) 

The important point is that |t) has been decomposed in terms of states it)} 
which do satisfy the time-translation property (Q) 

|t + r',±) = e^^°"'/^|t,±), (40) 
and transform under time reversal as 

n\t,±) = \-t,T), (41) 

so that the variable t, unlike r, could, in principle, be associated with physical time. 
Note that even though the states {|t,±)} constitute a complete set, 



they are not orthogonal, 



r+co 

/ dt\t,a){t,a\ = 1, (42) 

oo 



{t,a\t',a') = y / dE{t,a\E,P){E,P\t',a') = ^ dE e 

n: Jo h Jo 



OO 

iE{t~t')/h 



^6aa'ls{t-t')-P.P.—^\. (43) 



For this reason, the states it)} cannot be used to construct a self-adjoint oper- 



ator. 
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IV. MEAN ARRIVAL TIME 



One of the most controversial aspects of quantum mechanics is that concerning with 
the connection between the theoretical formulation and the corresponding measurement 
results. In its space-time representation, quantum mechanics becomes a continuous 
wave theory, whereas measuring devices usually deal with individual particles. The 
quantum formalism tells us how to obtain the probability distribution of the measure- 
ment results in terms of projections of the state vector onto appropriate subspaces of 
the Hilbert space. While in the standard interpretation it is commonly assumed that 
probability distributions refer to individual particles, their experimental verification 
requires an ensemble. Quantities defined in the ensemble may offer practical guidance 
not only in the interpretation of quantum measurement theory, but also in the search 
for the quantum counterpart of a classical physical variable. In this sense, the mean 
arrival time may be useful in looking for a probability distribution of arrival times. 

Consider a classical statistical ensemble of particles of mass m, directed along a well- 
defined direction, and characterized by the phase space distribution function f{x,p,t). 
The average time of arrival at a spatial point xq is given by 

/rcittj(xo,t) 

where J{xQ,t) represents the average current at xq, 

J{xo,t) = f{x,p,t) — 5(x — xo)dxdp, (45) 
J J m 

and plays the role of an unnormalized probability distribution of arrival times. It seems 
natural to make use of the correspondence principle in order to translate the expression 
for the classical average time of arrival, Eq. (|4^), to the quantum formalism. This 
can be accomplished by substituting J(xo,t) by the expectation value of the current 
operator 

J{X) = ^(P\X){X\ + \X){X\P). (46) 



2m 

Such a quantum definition for the average time of arrival has been widely used in 
recent years [|,18-2^. However, unlike its classical counterpart, even for wave packets 



directed along a well-defined spatial direction, the quantum probability current is not 
positive definite. For this reason, it cannot be considered as a probability distribution 
of individual arrival times, and the validity of the above expressions in the quantum 
context is questionable. In fact, strictly speaking, J{X) is an operator valued distribu- 



tion (the operator analog of a generalized function), and as pointed out by Wigner |25], 
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there exist quantum quantities, such as J(^), whose expectation values do not corre- 
spond to averages of individual measurements (eigenvalues) , but represent a measurable 
property of the ensemble whole. 

In spite of the general inadequacy of J{X) to describe the probability distribution 
of arrival times, when quantum backflow contributions become negligible the quantum 
current becomes positive and admits a probability interpretation |18|. Such a situation 
occurs, at large times, for freely moving packets containing only positive momenta, and 
it also occurs under the standard asymptotic conditions of scattering theory. In fact, Eq. 
(^) can be operationally justified in the quantum case by using a perfect absorber, i.e., 
a complex potential that absorbs the incoming wave in an arbitrary small spatial region, 
without reflection or transmission . According to such an operational model, which 
simulates the detection of incoming particles by a destructive procedure, the average 
time given by (^) coincides with the average time of absorption (detection) within any 
desirable precision. Thus, any properly defined arrival time probability distribution 
should be compatible with Eq. (^4|). 



V. PROBABILITY DISTRIBUTION OF ARRIVAL TIMES. 

FREE PARTICLE 

We shall restrict ourselves to the case of a free particle moving along a well-defined 
direction towards a detector situated at the point x = X. Specifically, we assume that 
the ingoing asymptote of the actual state of the particle corresponds, in the position 
representation, to a wave packet which is either a linear superposition of positive plane 
waves or a linear superposition of negative plane waves, 

IV'i.in) =©(±^)IV'±,in>- (47) 

Under these circumstances, the in asymptote becomes indistinguishable from the actual 
state \'ip±{t = 0)), so that we shall not discriminate between them from now on. Note 
that as a consequence of the commutation between the time-evolution operator e~*^o*/^ 
and the projectors G(ibP), at any time t it also holds 

|V^±(i))^e(±p)|V'±(t)). (48) 

As stated in the previous section, the mean arrival time at a point X is given by 

. . _ t^drr{i:^{r)\J{Xm{T)) 
/_rdr(V^±(r)|J(X)|V'±(r)) ' 
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where J{X) is the current operator in the Schrodinger picture, given in Eq. (^6|), and 
{iP±{t)\J{X)\'iP±{t)), is the probabihty current at the instant of time t = r, in the 
(Schrodinger) state \iP±{t)). 
In the free case, we have 

-oo 

dr(V'±(r)|J(X)|V±(r)) =±1, (50) 

so that, the mean arrival time {tx)± can be expressed as the following expectation 
value 

{tx)±^±{i^±\J±{XM±), (51) 
where \'4'±) denotes the state of the particle in the Heisenberg picture, i.e., 

|V;±)=e^^o-/^V±(r)) = |V±(0)), (52) 
and we have introduced the operator 

/+°° ^ ^ 

dTT@{±P)JH{X,T)e{±P), (53) 
-oo 

where Jii{X,t) is the Heisenberg current operator, 

Jh(X, t) = e'^°^/^ J{X) e~'^°^l^. (54) 
For later convenience use has been made in the above equations of the identity |V'±(0)) = 

e(±p)|v^±(o)). 

Inserting twice the resolution of unity, Eq. (p!4[), and using 

Q{±P)\E,a) = 5a,±\E,a), (55) 

J'±{X) takes the form 

dTT dE dE'e'^^-^^^/^\E,±){E,±\J{X)\E',±){E',±\. (56) 

-oo Jo Jo 

Substituting the expression (^) for J{X), using Eq. (|T^), and taking into account that 



according to Eq. (16) 



{X\E, ±) = h~^/^{m/2E)^/^ ^±iV2^x/n^ (57) 
the matrix element in the integrand of (p6|) can be rewritten in the form 
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{E, ±\J{X)\E', ±) = ±^|(|)'^'+ (f )'^'} e^<v/^-^^)^/^ (58) 

After insertion of Eq. (|5^ ) into Eq. (^), the operator J^iX) reads 

J±{X) = e-'^^/^ J±{X = 0) e+'^^/^, (59) 
where we have again taken advantage of Eq. ( p!3| ) to write 

and J±(^ = 0), which is the operator involved in the determination of the mean arrival 
time at the point X = 0, is given by 



. ./>.p.p.i^{(|)-'% (I)-} 

(61) 

In order to guarantee that the integrand is well behaved over the whole interval 
of integration we shall restrict ourselves to physical states satisfying the boundary 
conditions 

lim E^/"^ {E,±\ip±) =0, limS~i/^(^,±|'0±) = 0, (62) 

which, in the more familiar momentum representation, take the form 

lirn {p\Tp±) = 0, lini p'^ {p\ip±) = 0, (63) 

p— >±co p-+0 

Put another way, we shall restrict ourselves to normalizable wave packets, which are 
superposition of either positive or negative plane waves, and vanishing faster than p as 
p approaches zero. 

The integral in the r variable, in Eq. (|6l|), can be readily performed to obtain 

A(0) = 4f.EpE'(|^*(E-E')){(|)"%(|)"]|E.±>(E'.±|. (04) 

To proceed further it is convenient to consider the matrix elements of i7±(0) between 
arbitrary states |$), l^*) satisfying the boundary conditions (|6^). Using the derivative 
of the Dirac delta in the integrand of (^) to perform one of the two energy integrals, 
one arrives at 
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=±ifij^ dE{E,±\^)—{^\E,±). (65) 



On the other hand, using the resolution of the unity, Eq. (32), and taking into account 
that 

±in-^{T\E,±)=T{T\E,±), (66) 

one can obtain a useful alternative expression for the energy derivative in the integrand 
of Eq. (|65|), in terms of the self-adjoint "time" operator defined in Eq. (^). Indeed, 

±ih — {^\E,±) = ±ihJ_^dT{'^\T)—{T\E,±) 

/ + 00 
dTT{^\T){T\E,±) = {<S>\T\E,±). (67) 
-oo 

Therefore, 

POO 

($|J±(0)|^) = / dE{^\f\E,±){E,±\^). (68) 
Jo 

Taking |$) = |^) = 6+*^-^/^ | ^± ) , we have 

{tx)± ^ ±(V'±|e-^^^/^J±(0)e+*^^/^|V±) 

roo 

= ± dE{'ilj±\e-'^^/^f\E,±){E,±\e+'^^/^\^l;±). (69) 

^0 

Using the identity 

\E,±){E,±\ ^ Q{±P) ^ \E,a){E,a\, (70) 

as well as Eq. (0), and taking into account that e(±P) e+*^^/^|V'±) = e+^^^/^\^±), 
we finally find 

{tx)± = ±(^±|e-^^^/^te+*^^/^|^±). (71) 

Accordingly, the self-adjoint operator involved in the determination of the mean 
arrival time at an arbitrary point X, is given by the spatial translation of the operator 
T previously defined. 
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t(X) = e-*^^/'*re+^^^/^ (72) 
and its spectral resolution reads 

/ + 00 
dTT\T;X){T-Xl (73) 

where 

r'+OO 



r;. 



Since the states {|r;X)} are generated from the complete and orthonormal set 
{|t)} via a unitary transformation, they also constitute, for a given X, a complete and 
orthonormal set. 

Introducing the complete but nonorthogonal set of shifted states {\t, it; X)}, defined 
as the spatial translation of the set it)}, 

/+00 , 
dEe'(^'^'^''y''\E,±), (75) 
-00 

the states \t;X) can be decomposed as a superposition of negative- and positive- 
momentum contributions, in the form 

\t;X) = \t = -T,-;X) + \t = +T,+-X). (76) 



Inserting now the spectral resolution of T{X), Eq. ([73D, into Eq. ( [TlD one can 
express the mean arrival time at the spatial position X, in the form 

/ + 00 
(±T) KT;X|V±>|^ (77) 

so that, 

/ + 00 
drTK^ = +r,+;X|V+>|^ (78) 
-oo 



+00 



{tx)- = I dr i-r) \{t = -T,--X\i;^)\\ (79) 
On the other hand, taking into account the resolution of the unity, 

/-(-oo 
dt \t,a;X){t,a;X\ = 1, (80) 
-oo 



a=±'' 
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we have 



/+0O 
dt\{t,±-X\xlj±)\\ (81) 
-oo 

which, for a free particle, coincides with the total probability of arriving at the point 
X at any instant. 

Therefore, the quantities |(t;X|V'±)P enter the above equations as a probability 
density, and lead to an expression for the mean arrival time having the correct semiclas- 
sical limit in terms of the probability current. However, unlike the latter, it is definite 
positive. Accordingly, for a free particle in the Heisenberg state one can interpret 
consistently (t;X|'0+) = (t = +r, as the probability amplitude of arriving at 
the spatial point X from the left, at the instant t = t. Similarly, for a free particle in 
the Heisenberg state \tp-), the scalar product (— r;X|'0_) = {t = +t, —; X\ip-) can be 
interpreted as the probability amplitude of arriving at X from the right, at the instant 
t = T. 



VI. PROBABILITY DISTRIBUTION OF ARRIVAL TIMES. 
POTENTIAL BARRIER 

Consider the passage of particles incident from the left over a one-dimensional 
potential barrier V{x). As usual, we assume that far away from the scattering center, 
V{x) vanishes sufficiently fast as to guarantee the validity of the standard scattering 
theory formalism. Under the conditions we are interested in, the ingoing asymptote, 
I'^in), of the actual state of the particle satisfies 

IV'in) = e(P) l^in). (82) 

The M0ller operators which play a central role in scattering theory, are defined 

as 

n±= lim e^^*/^e-'^o*/^, (83) 

where Ho = P'^ /2m, and H = Hq + V{X) is the Hamiltonian governing the dynamical 

evolution of the system. These operators have the importance that they map the 
asymptotic states onto the corresponding scattering states. Specifically, the actual 
state of the particle, \'ip{t = 0)), is related to its in and out asymptotes, iV'in) and iV'out)) 
by means of 



16 



m = 0)) =n+\tl;i^) = n_\ij,^i). (84) 
Making use of the intertwining relations for the M0ller operators ||2^ , 

nlHQ± = Ho, (85) 

the mean arrival time at a spatial point X, on the right of the barrier and asymptotically 
far from the interaction region is given by 



/_+-drr(^(r)|J(X)|^(r) ) 
J^::dT{^p{T)\J{X)\^l;{T)) 

where now we have 



(V(r)|J(X)|V(r)) = (V-inle^^^^/^ J^^ J(X)f]+e-*^o^/^Vin>. (87) 

Inserting twice the resolution of unity, Eq. (0), and taking advantage of Eq. (|8^ ) 
to write (pltpm) = @{p) {p\ipin), one obtains 



(V'(r)|J(X)|V'(T)) = / dp' dpe^^p-/^-*^W?^(Vin|p')(p'|j^U(^)^^+|p)(p|V'in), 
Jo Jo 

(88) 

where Ep = /2m. 

The state \p+) = O+lp), which is the solution of the Lippmann-Schwinger equation 
corresponding to an ingoing plane wave \p), satisfies the eigenvalue equation H \p+) = 
Ep \p+) with the boundary conditions 

X — > —00 : {x\p+) ~ {x\p) + R{p) (x| —p), (89) 



X — > +00 : ~ r(p) (x|p), (90) 

R{p) and T{p) being the reflection and transmission coefficients, respectively. Thus, 
for spatial points X on the right and asymptotically far from the interaction center one 
has 

{p'+\J{X)\p+) = T*{p')T{p){p'\J{X)\p), (91) 
where use has been made of Eqs. (p6D and (pO|). Substituting in (38) we obtain 
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(V(T)|J(X)|V^(r)) =/ dp' dpr*(p')(V'in|p')(p'|e^^o^/^J(X)e-^^o-/^|p)(p|V'i„)T(p). 
Jo Jo 

(92) 

Taking ( |92|) into account, the r-integral in the denominator of (^) can be readily 
carried out to obtain a Dirac's delta, 

5{E.p, -E,) = - 5{p' -P) + ^ 6{p' + p), (93) 
\P\ \P\ 

and this can in turn be used to finally obtain 

dT{i,{T)\J{XmT))= dp\T{p)\^\{p\i;i^)\\ (94) 



which is nothing but the transmittance. Defining the (unnormalized) freely evolving 
transmitted state, \iptr), as 



= / dpT{p){p\ij,^)\p), (95) 
Jo 

and using Eq. (|9^), the mean arrival time (|8^) takes the form 

(tx) = I , , (V'trl / dTre(P)Jl(X,T)G(P)|Vtr.), (96) 

where we have taken advantage of Eq. (^2|) , and Ji(X, r) denotes the current operator 
in the interaction picture, 

Ji{X, r) = e^^o^/^ J{X) e~'^°^/^. (97) 



A comparison between Eqs. ( p7\ ) and ( |54D shows that the current operator in the 
interaction picture coincides with the free current operator in the Heisenberg picture. 
This is the important point which allows us to rewrite {tx) in terms of the freely evolving 
operator J'^(X) of Eq. ( |53D and, consequently exploit the formalism developed in the 
previous section for the free case, to finally obtain 

^^^^ ^ iiJtrlMXmr) ^ {i^tr\f{X)\iJtr) ^ ^gg^ 
ii'trli'tr) {tptr\tptr) 



Inserting the expression for the self-adjoint operator T{X), given in Eqs. (|7^) and (74), 
the mean arrival time at a spatial point X behind the barrier and asymptotically far 
from the interaction center, takes the suggestive form 
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1 r+oo 

{tx) = , , I , , / dTT\{T;X\i;tr)\^. (99) 



Furthermore, taking into account that, for a given X, the states {|t;X)} constitute a 
complete and orthonormal set, we find that the integral 



oo 



dT\{T;X\i;tr)\' = {i;tr\i^tr)= dp\T{p)\'\{p\i;,^)\', (100) 

D Jo 

coincides with the transmittance, which is nothing but the total probability of arriving 
at an asymptotic point behind the barrier. Therefore, we can consistently interpret 
{T;X\iptr) = {t = T,+; Xliptr) as the (unnormalized) probability amplitude of arriving 
at the asymptotic point X, behind the barrier, at the instant t = t. 

The above results can be expressed in terms of the ingoing asymptote iV^in); by using 
the scattering operator 5 = qI_ which relates the in and out asymptotes, IV'out) = 
5 \ ipin)- Indeed, it is shown in the Appendix that the freely evolving transmitted state 
\tptr) can be written as 

|Vtr.) = e(p)5|Vi„), (101) 

so that the (unnormalized) probability density of arriving at an asymptotic point X, 
behind the barrier, at the instant t = t reads 

|(r; X\ e{P) S |Vin)|' = \{t = r,+;X\S |V'in)|'. (102) 

Finally, it should be noted that when the wave packet corresponding to the actual 
scattering state at i = (which is assumed to be a linear superposition of only positive 
plane waves) does not overlap with the potential barrier, then it becomes physically 
indistinguishable from the asymptotic ingoing wave packet and the above equations 
hold true with the substitution \ipin) \tp{t = 0)). 



VII. TIME-ENERGY UNCERTAINTY RELATION 

Giving a precise meaning to the well-known time-energy uncertainty relation seems 
to be a reasonable requirement for any quantum formulation of the time-of-arrival 
concept. 

As already said, the commutation relation [H, T] = ifi automatically leads to the 
uncertainty relation (|3|). Although, it has not been possible to develop a quantum 
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formulation of the arrival-time concept based on such a commutation relation, there 
is still room for a time-energy uncertainty relation, because even though the existence 
of a commutation relation is a sufficient condition for the existence of an uncertainty 
relation, it is by no means a necessary condition. 

It should be noted that the self-adjoint operator T(X) defined by Eqs. (^)-([7^) is 
conjugate to the operator Id. Specifically, 

[n, f{X)] = e-*^^/^ [TY, t] e+'^^'^ = in. (103) 

Introducing a probability amplitude for the time of arrival, in terms of the eigenval- 
ues and eigenstates of a self-adjoint operator satisfying the above commutation relation 
has as an important consequence the existence of a time-energy uncertainty relation. 
To see this, let us consider the problem studied in the previous section, namely, the 
arrival time of a quantum particle at a detector located behind a potential barrier and 
asymptotically far from the interaction center. [The free case is nothing but a particular 
case of the latter corresponding to T{p) — > 1, which implies \'4^tr) — *■ iV'in)-] 

Because of Eq. ( |103D , it automatically holds 

AW ATx > n/2, (104) 

where A'H and ATx are the root-mean-square deviations of the corresponding observ- 
ables,^i.e., {^Uf = (H^) - {Hf, and (ATxf = {f^{X)) - {f{X)f, with (A) = 
{iptr\A.\il>tr) / ii^trlfptr) ■ However, from Eqs. ( p3| ) and ( |95| ) it follows that 

n\^Ptr) = Hol^Ptr), (105) 

and hence ATI coincides with the statistical spread of the energy of particles arriving 
at the detector, 

(Anf = {Hi) - {Hof ^ {AEf. (106) 
On the other hand, according to Eqs. (|73|) and (|98| ) we have 

{ATxf = {{f{X)-{tx)) ) = — — ^ / dr{r-{tx)?\{T-X\i„r)\^^{Atx)\ 

(107) 

and since \{t; X\iptr)\'^ / {iptr\'4'tr) is the probability distribution of the arrival time of 
particles at the detector, the above equation shows that ATx is nothing but the cor- 
responding statistical deviation, Atx- Therefore, the statistical spreads of the energy, 
AE, and time of arrival, Atx, of particles reaching the detector satisfy the time-energy 
uncertainty relation 
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AE Atx > fi/2. (108) 



VI. CONCLUSION 

Despite its fundamental nature, the quantum formulation of the time-of-arrival 
concept is a problem which remains open nowadays. This question has the additional 
interest that probability distributions of arrival times are, in principle, experimentally 
accessible via the time-of-flight technique. Moreover, a quantum formulation of such 
a problem may provide a useful tool for a better understanding of the tunneling time 
problem as well as its possible technological applications. 

The main difficulty for defining a quantum time operator lies in the nonexistence, in 
general, of a self-adjoint operator conjugate to the Hamiltonian, a problem which can 
be always traced back to the semibounded nature of the energy spectrum. In turn, the 
lack of a self-adjoint time operator implies the lack of a properly and unambiguously 
defined probability distribution of arrival times. 

Although it has been shown that under certain circumstances of physical inter- 
est, the probability current becomes positive and admits a proper interpretation as an 
unnormalized probability distribution of the time of arrival []l^20|, it cannot be con- 



sidered as a fully satisfactory solution of the problem, for it is not positive definite as 
it should be. 

In searching for a probability distribution defined through a quantum time operator 
one has to circumvent the problem stated above. There are two possibilities. If one de- 
cides that any proper time operator must be strictly conjugate to the Hamiltonian, then 
one has to renounce to find a self-adjoint operator. (Even though such a property is a 
hallmark of any observable in the standard quantum formalism, it is not strictly nec- 
essary for a consistent formulation of probability distributions of measurements results 



|27,2^.) If, on the contrary, one imposes self-adjointness as a desirable requirement 
for any observable, then one necessarily has to abandon the requirement that such an 
operator be conjugate to the Hamiltonian. In the present paper we have adopted the 
latter via. We have explicitly constructed a self-adjoint operator Ti. with dimensions of 
energy, and a non-bounded spectrum. Such an operator is essentially the energy of the 
particle with the sign of its momentum. The non-bounded character of its spectrum 
enables us to introduce a self-adjoint operator with dimensions of time, T, by demand- 
ing it to be conjugate to Ti.. Since the latter is essentially the Hamiltonian, except 
sign, one expects the self-adjoint operator T so defined to be physically meaningful 
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and relevant to the arrival time problem. Indeed, we have shown that it is possible to 
define consistently a probability distribution of arrival times at a spatial point, in terms 
of the eigenvalues and eigenstates of such an operator. This probability distribution, 
which is a function of the initial state of the system, does not depend on the particular 
design of the measuring device, and has the additional desirable consequence that it 
leads to a precisely defined time-energy uncertainty relation. 



APPENDIX: TRANSMITTED STATE AS A FUNCTION OF THE 

INGOING ASYMPTOTE 

In this appendix we show that the freely evolving transmitted state \iptr) can be 
written in terms of the ingoing asymptote \^p■m) in the form [Eq. ( |101 )] 



\lPtr) = e(P) IVout) = e(P) S IVin), (109) 

where S denotes the scattering operator, relating the ingoing and outgoing asymptotes, 

/+00 r+oo 
dp dp\p'){p'\S\p){p\^Pin). (110) 
-oo J — oo 

Taking into account that the matrix elements of the scattering operator can be 



written in terms of the on-the- energy- shell T matrix as [^| 

{P\ S \p) = 6{p - p') - 27ri 5{Ep - E^,) {p'\ J{Ep + iO) \p), (111) 
equation (|110|) reads 



+ 



IV'out) = S \t,n) = 1^ dp (l - 2Trim/p {p\ J{Ep + iO) \p)) {p\i^i^) \p) 

dp [-2TTim/p{-p\JiEp + iO)\p)j (plV^in) \-p). (112) 



On the other hand, from the Lippmann-Schwinger equation for \p+) it follows that 
the wave function {x\p+) can be written as 

{x\p+) = {x\p) + J dx'{x\ {Ep + iO- Hoy^\x'){x'\ J{Ep + iO) \p). (113) 
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Substituting in the above equation the expression for the Green's function (which can 
be easily obtained by inserting the resolution of unity in terms of momentum eigenstates 
and evaluating the resulting integral by contour integration in the complex plane), 

(x| {E„ + iO- Hoy^\x') = -^e^lPll^-^'IM^ (114) 

h\p\ 

one obtains, for p > and x — > +00, 

{x\p+) ~ {x\p) — 2-111171 /p {x\p) J dx' {p\x'){x'\T{Ep + iO) \p) = 

(1 - 2TTim/p {p\ J{Ep + iO) {x\p). (115) 

A comparison with Eq. ( |90| ) yields 

T{p) = (1 - 2Trim/p {p\ J{Ep + iO) . (116) 
Similarly, for p > and x — > — cxd, Eq. ( \113\ ) leads to 



{x\p+) ~ {x\p) — 2iTim/p {x \—p) J dx' {—p\x'){x'\ T{Ep + iO) \p) = 

{x\p) + (^-2TTim/p {-p\ J{Ep + iO) {x\p-), (117) 

and comparing with Eq. (p9D we find 



R{p) = i^-2mm/p {-p\ J{Ep + iO) |p) j . (118) 
Substituting Eqs. ( 116p and ( |118| ) in ( |112p we finally arrive at 



/•oo rco 

\^ont) = S\t.n)= dpT{p){p\i,,^)\p)+ dpR{p){p\tm)\-p), (119) 
Jo Jo 



from which it follows Eq. (109). 
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